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Abstract 

Solar-like oscillations are stochastically excited by tur- 
bulent convection. In this work we investigate changes in 
the acoustic oscillation power spectrum of solar-type stars 
by varying the treatment of convection in the equilib- 
rium structure and the properties of the stochastic excita- 
tion model. We consider different stellar models co mputed 
with the stand ard mixing-length description by Bohm 
Vitense (1958) and with a generalized formulation of the 
mixing-length approach by Cough (1976, 1977). We cal- 
culate the acoustic power generated by the turbulent con- 
vection which is injected stochastically into the acoustic 
pulsation modes. 

Large differences in the oscillation powers are obtained 
depending on the choice of the assumed convection for- 
mulation. We show that the high-quality data Eddington 
will provide, will allow us to distinguish between theoret- 
ical predictions of acoustic power spectra obtained with 
different convection models. 



1. Introduction 

The amplitudes of solar-like oscillations are defined by the 
balance between excitation and damping. The oscillations 
are excited stochastically by the acoustic noise generated 
by the turbulent motion of the convective elements. In 
Samadi & Goupil (2001, Paper I hereafter) a theoretical 
formulation for the acoustic power injected into solar-like 
oscillations is proposed and which supplements previous 



theories. We refer the reader to Samadi (2001) for a de- 
tailed summary and discussion on some recent unsolved 
problems. 

The excitation process depends on the assumed turbu- 
lence spectrum, as discussed, for example, by ^amadi et al 



(2001b) and 3amadi (2001). It also depends crucially on 



the convection model to compute the stratification of the 
convectively unstable layers in the equilibrium model. The 
amount of energy injected into the oscillations depends 
strongly on the velocity of the convective elements. 

The main goal of this work is to asses changes in the 
oscillation power spectrum due to modifications of the 
convection treatment in the equilibrium model. We con- 
sider two different formulati ons: the classical desc ription 
of mixing-length theory by Bohm- Vitense (1958) and a 



nonlocal generalization of the mixing-length formulation 
by Gough (1976, 1977). Additionally we study the depen- 
dence of the oscillation power on the assumed turbulence 
spectrum in the excitation model for both convection for- 
mulations. 

The theoretical formulation in Paper I involves two free 
parameters. These parameters are calibrated such as to 
reproduce for a solar model the observed acoustic power 
spectrum (Section 2). In Section 3 we compute oscillation 
power spectra for several stellar models using the two con- 
vection formulations. 

We conclude that Eddington's performance will allow 
us to distinguish between the two treatments of convec- 
tion and consequently that Eddington will provide further 
constraints on convection models. 

2. Power injected into solar-like oscillations 

2.1. Theory of the stochastic excitation 

The acoustic power P injected into the oscillations is de- 
fined {e.g., Goldrcich ct al. 1994) in terms of the mode 



damping rate rj, the oscillation mean-square amplitude 
(v4^), the mode inertia / and oscillation frequency uj as : 

Piu) ^ rj (A^) I . (1) 

The mean-square amplitude is defined by both the ex- 
citation by turbulent convection and by the damping pro- 
cess. It can be written as 

X {§j^,{L0,m)+3i\m):¥^^„m)§s{L^,m)} , (2) 

where is the radial displacement eigenfunction, p the 
density, A is the mixing length, w the vertical rms velocity 
of the convective elements; 3^{£,r,m) is a function which 
includes the second derivative of ^r, and 31^ (to) describes 
the ratio of the excitation by the entropy fluctuations to 
that by the Reynolds fluctuations. The source functions 
§^{LL),m) and §s(a;,TO) describe the contributions from 
the Reynolds and entropy fluctuations, respectively, aris- 
ing from the smaller scales of the turbulent cascade. De- 
tailed expressions for {A'^'j, §^{uj,m), §s{oj,m), 51^ and 
3^^ are given in Paper I. 

The source functions include the turbulent kinetic en- 
ergy spectrum E(k), and the turbulent spectrum of the 
entropy fluctuations Es{k) which can be related to E{k) 
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by a simple expression {e.g., Samadi et al. 2001a, Paper II 
hereafter). Both source functions SR(a;,m) and §g(ti;,m) 
are integrated first over all eddy wavenumbers k, followed 
by an integration over the stellar mass M to obtain the 
acoustic power P (see Eq]|). 

In the present work, p and w are obtained from the two 
equilibrium models, computed with the aforementioned 
two convection formulations. The eigenfunctions and 
eigenfrequencies lo are obtained from two different oscilla- 
tion codes and the fc-dependence of E(k) is inferred from 
different observations of the solar granulation and from 
theory. 

2.2. The equilibrium models 

We consider two sets of stellar models: the first set con- 
sists of complete models obtained with the CESAM evo- 
lutionary code in which the convective heat flux is com- 
puted according to the classical mixing-length theory by 
Bohm-Vitense (1958, C-MLT hereafter). The momentum 
fiux (sometimes referred to as turbulent pressure) is ne- 
glected in this set of models. The eigenfunctions are ob- 
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Figure 1. Kinetic turbulent spectra versus scaled wavenum- 
ber fc/fco- 

cade, the spectra exhibi t different scaling laws. These spec- 



tra were considered by pamadi et al. (2001b)| to compute 



tained from the adiabatic FILOU pulsation code by Tran 



acoustic power spectra for various solar-type stars and are 
also considered in this paper. 



Minh ^. Lonu n C)Q5^|. The detailed input phvsics used in From ob serv ations of the solar granulation [Espagnet 



this set of models is described in Paper IT. 



et al. (1993) and Nesis et al. (1993) suggest various values 



The second set consists of envelope models computed 



in the manner of Balmforth (1992a) and Houdek et al, 
(1999). In these models convection 



is 



treated with the 
nonlocal mixing-length formulation by Gough (1976, G- 
MLT hereafter), which consistently includes the momen- 
tum flux (i.e., the rr-component of the Reynolds stress 
tensor). The eigenfunctions are obtained from the nona- 
diabatic pulsation code by Balmforth (1992a), which in- 



cludes both the Lagrangian perturbations of the heat and 



a H. 



momentum fluxes in the manner of Gough (1977) 
For both model sets the mixing length A = 
where Hp is the local pressure scale height and a the 
mixing-length parameter, is calibrated first to a solar model 
to obtain the helioseismically determined depth of the con- 



vection zone of 0.287 solar radii (Christensen-Dalsgaard 
et al. 1991D. 



2.3. Models for stellar turbulence 

Several turbulent spectra E{k) were discussed in Paper II: 
the "Nesis Kolmogorov Spectrum" (NKS hereafter) and 
the "Raised Kolmogorov Spectrum" (RKS hereafter) were 
suggested from observations of the solar granulation by 



Espagnet et al. (1993) and Nesis et al. (1993), Here we 
also consider the "Broad Kolmogorov Spectrum" (BKS 
hereafter) by Musielak et al. (1994). These spectra are de- 
picted in Figure For wavenumbers k > ko, where fcg 
is the smallest wavenumber of the classical Kolmogorov 
spectrum, all spectra follow the classical Kolmogorov scal- 
ing law, k~^^^. Only in the low- wavenumber range, k < ko, 
where kinetic energy is injected into the turbulent cas- 



tor fco- Because of this ambiguity in the choice of ko, we 
relate ko to the mixing length A by /cq = 27r / (/3A), where 
/? is a free parameter of order unity (Paper I). Another un- 
certainty in the formulation of the excitation model it the 
eddy correlation time scale. Because of our still poor un- 
derstanding of turbulent convection the eddy correlation 
time scale is not well defined and therefore needs to be 
scaled, leading to an additional free parameter A of order 
unity {e.g., Balmforth 1992b). As it was demonstrated in 
Paper II the oscillation power computed for the Sun de- 
pends crucially on the values of the free parameters A and 
p. 

2.4. Calibration of the free parameters 

The mean square surface velocity Vs is related to the acous- 
tic power P{lu) and the mode damping rate rj by the ex- 
pression : 

v^,=eirs)P{u^)/2vI, (3) 
where rg is the radius at which the oscillations are ob- 
served. 

The acoustic power P is computed for the two cali- 
brated solar models using the aforementioned convection 
formulations. For the velocity estimates Vs, Eq.^, we as 



sume the observed linewidths (damping rates) by Lib 



brecht (1988). The free parameters f3 and A are calibrated 
for all turbulent spectra (see Section 2.3) and for both 
solar models to fit the estimated velocities Vs to the ob- 
servations by p^^ibbrccht (1988) . 



Figure || displays the results for the solar model com- 
puted with G-MLT. The corresponding calibrated values 
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Figure 2. Computed surface velocities Vs (Eq]^ assuming 
various turbulent spectra (see Fig.^ and G-MLT for com- 
puting the turbulent fluxes. The free parameters A and (3 
are calibrated to Libbrecht's (1988) velocity and linewidth 
(damping rate) measurements. 

Table 1. Calibrated values of the parameters [3 and A for 
the solar model computed with C-MLT. 



spectrum 


p\ 


P 


A 


RKS 


0.6 


1.96 


0.31 


BKS 


1.6 


4.06 


0.39 


NKS 


2.6 


3.11 


0.83 



of (3 and A are listed in Table |i| (see Paper II for the cali- 
bration results of the models computed with the C-MLT). 

At low frequencies the velocity Vs computed with G- 
MLT is in better agreement with the data than the results 
published in Paper II obtained with the C-MLT. At high 
frequencies the differences in Vs obtained with the various 
turbulent spectra are of similar small magnitude than the 
results of Paper II which assumed the C-MLT. Best agree- 
ment between theory and measurements is obtained with 
the NKS for both convection formulations. 

3. Scanning the HR diagram 

3.1. Stellar models 

We consider various solar-type stars with masses between 
1 Mq and 2 Mq in the vicinity of the main sequence. 
The model parameters are listed in Table S, which cor- 



respond to the models considered previously by Samadi 
et al. (2001b)| . As for the solar models in Section 2 we 
consider two sets of stellar models: the first set, computed 
with the C-MLT, is adopted from ^amadi et al. (2001b) 



The second set, computed with G-MLT, assumes the cali- 
brated mixing length of the solar model discussed in Sec- 
tion 2. The acoustical cut-off frequencies are very similar 
between the two sets of models and are displayed in Ta- 
ble for the first set. 



Models 


L 




M 


Age 








m 


[Mq] 


[Gyr] 


[mH 


A 


12.1 


6350 


1.68 


1.79 


1.0 


B 


9.0 


6050 


1.44 


3.05 


1.0 


C 


6.6 


6400 


1.46 


2.40 


1.5 


D 


3.7 


5740 


1.08 


7.33 


1.5 


E 


3.5 


6120 


1.25 


4.10 


2.3 


F 


2.6 


6420 


1.25 


1.76 


3.6 



Table 2. Model parameters of solar-type stars. The model 
age and acoustical cut-off frequency Vc «re listed for the 
models obtained with the CESAM evolutionary code which 
assumes the C-MLT for computing the convective heat 
flux. 



3.2. Dependence on convection models 

For all stellar models in both sets we assume the NKS and 
the calibrated values of (3 and A quoted in ^amadi et al 
(2001b)| and in Table 0. 



The position in the HR diagram of all stellar models 
and a qualitative overview of their acoustic power spectra 
P are depicted in Fig. |[ Detailed results of P are shown 
in Fig. ||. 

At high frequencies the differences in P between mod- 
els computed with the C-MLT and G-MLT increase with 
increasing effective temperature Tos and luminosity L. As 
discussed in Houdek (1996) , the nonlocal formulation (G- 
MLT) predicts smaller temperature gradients in the upper 
superadiabatic region relative to the C-MLT. This means 
that convection is more efficient in models computed with 
G-MLT which leads to a different profile (depth-dependence) 
of the superadiabatic temperature gradient between the 
C-MLT and the G-MLT. The differences in the superadi- 
abatic temperature gradient between the two model sets 
increase with L and Toff. These results are illustrated in 
Fig. 1^ which shows the superadiabatic temperature gra- 
dient V — Varf versus i?* — r, with i?* being the radius 
of the star. With increasing L or Tefi, the maximum of 
V — Vad is shifted more rapidly to deeper layers for the 
G-MLT models than for the C-MLT models. This leads 
to progressively larger differences in the convective veloc- 
ities w and in the shape of the eigenfunctions between the 
two sets of models. These differences in w (note that P 
depends crucially on w, see Eq.^) and in the shape of the 
eigenfunctions are particular large in the superficial layers 
and result in a larger amount of acoustic power injected 
into high frequency modes for models computed with the 
C-MLT. 

There is an additional age effect: it can be seen from 
comparing the results of the models E and F (see Fig. ||) . 
Model E has the same mass than model F but is older. 
The increase of the maximum acoustic power with age is 
found to be larger for models computed with the C-MLT 
than for the G-MLT models. 
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Figures. Positions of models in the HR diagram (filled circles). For each model of the two sets, the computed oscillation 
power spectra assume the NKS and are illustrated in an associated box : the solid curves correspond to models computed 
with G-MLT and the dashed curves display the results obtained with the C-MLT. The plot range of the abscissa starts 
from zero and extends to the acoustic cut off frequency, indicated for each model. The power spectra are normalized 
with their maximum values. 



3.3. Dependence on turbulence models 



3.4. Discussion 



We compute P for all the kinetic turbulent spectra dis- 
cussed in Section 2.3. The results are shown in Fig.^ for 



One of Eddington's tasks will be the continuous observa- 
tion of stellar luminosity oscillations over a time period 



of 30 days with a frequency accuracy of 0.3 /iHz (Favata 



the C-MLI' (top panels), the dependence of P on the as 
sumed turbulence spectra is more pronounced for the hot- 
ter model C than for model E. This is not observed for 
the model set computed with G-MLT (bottom panels). 

Models computed with the C-MLT exhibit a more pro- 
nounced decrease of the depth of the surface convection 
zone with increasing Tcff. This leads to a smaller extend 
of the excitation region for hotter models compared to 
the models computed with G-MLT. As discussed in more 
detail by Samadi et al. (2001b, c) a shallower excitation re- 
gion results in a more pronounced frequency-dependence 
of P on the assumed turbulent spectra. 

Model C, obtained with G-MLT, exhibit a deeper ex- 
citation region and consequently the dependence of P on 
the assumed turbulent spectra is smaller. This property 
can be explained by the larger efficacy with which G-MLT 
transports the convective heat flux. 



t |ie models C and n.. for the set oi models computed with 2000| ). Furthermore the large telescope of Eddington 



(approximately ten times larger than that of the COROT 
mission) will lead to a noise level of ~ 0.2 ppm for stars 
with a magnitude = 6 and for an observing period of 
30 days. Fo r comparison, COROT will reach a n oise level 
of 0.7 ppm (lAuvergne fc the COROT Team 2000|) for stars 
with a magnitude mi, ~ 6 and for a continuous observing 
period of 5 days. 

Will the high-quality observations of the Eddington 
mission be accurate enough to allow us to distinguish be- 
tween the changes in the predicted oscillation powers P of 
Section 3.2 and Section 3.3 ? In order to answer this ques- 
tion we estimate the expected accuracy of Eddington's 
measurements of P. 

Kjeldsen & Bedding (1995) suggested a very simplified 
scaling law between oscillation amplitude ratios {SL/L) / Vg 
and Toff: 

6L/L <x Vs tJ'^ , (4) 

assuming adiabatic oscillations in a purely radiative star. 
It suggests that the amplitude ratios scale inverse propor- 
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Figure 4^. Oscillation power spectra computed for the models C and E assuming the NKS (solid curves), the RKS (dashed 
curves) and the BKS (dot-dashed curves). The top panel displays the results for models computed with the C-MLT and 
the bottom panel illustrates the results for models computed with G-MLT. Vertical error bars AP are obtained from 
Eq.^) assuming an accuracy for rj/ 2t: (resp. 5L/L) of 0.3 fiHz (resp. 0.2 ppm). 



tionally to the effective temperature (but see also Houdelc 
et al., 1999, Fig. 14, who found that the amphtude ratios 
scale directly with TcfA. _ 

According to Eq.(|^) and Eq, 
P can be expressed as 
A77 



the relative error for 



AP A(5L 
= 2-^ — 

P SL 



V 



(5) 

where SL is obtained from Eqs.(^, ^). The damping rates 
77 are supplied fro m the nonadiabatic pulsation code of 
Balmforth (1992a) assuming equilibrium models computed 
with G-MLT. 

In Figure m and I the error bars AP are plotted accord- 
ing to Eq.(^) assuming a noise level of A{SL) ^ 0.2 ppm 
and A(7;/27r) ~ 0.3 /iHz. 

We conclude that for hotter stars with a magnitude 
ruy < 6, Eddington will provide data of sufficient accu- 
racy which will allow us to distinguish between the results 
obtained with the two considered convection formulations 
and consequently will provide further details on how to 
improve stellar convection models. 
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Figure 5. Oscillation power spectra computed for both sets of stellar models. The NKS turbulence model is assumed. 
The solid curves display the results for models computed with G-MLT and the dashed curves show the results for models 
computed with the C-MLT. 




Figure 6. Depth dependence of the superadiabatic temperature gradient V — Vad for the two sets of stellar models with 
surface radius i?* . The solid curves display the results for models computed with G-MLT and the dashed curves are 
the results assuming the C-MLT. 
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